On the Classical Model for Microwave Induced Escape from a Josephson Washboard 
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We revisit the interpretation of earlier low temperature experiments on Josephson junctions under 
the influence of applied microwaves. It was claimed that these experiments unambiguously estab- 
lished a quantum phenomenology with discrete levels in shallow wells of the washboard potential, 
and macroscopic quantum tunneling. We here apply the previously developed classical theory to a 
direct comparison with the original experimental observations, and we show that the experimental 
data can be accurately represented classically. Thus, our analysis questions the necessity of the 
earlier quantum mechanical interpretation. 
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Since early days, Josephson junctions in superconduct- 
ing circuits have been modeled very successfully with the 
so-called Resistively and Capacitively Shunted Junction 
(RCSJ) equivalent circuit 1]. In the situation where 
a junction is biased with a fixed current, that equiva- 
lent circuit yields an equation of motion for the junction 
phase that is analogous to a fictitious particle confined 
to a "washboard" potential consisting of repeated wells. 
In the presence of thermal noise, a Josephson junction 
biased just below its critical current may be excited out 
of its position at the bottom of a well and into a running 
state with an accompanying voltage. In 1981, a paper 
appeared f2| which reported experiments carried out on 
such a current biased Josephson junction, but at much 
lower temperatures (down to 3mK) than previously. It 
was found that at high temperatures the transition rate 
out of the metastable state was well described by ther- 
mal activation over the barrier formed by the rim of the 
well, but that as T — )■ 0, the transition rate was dom- 
inated by quantum mechanical tunneling of the macro- 
scopic junction phase (MQT). The transition from clas- 
sical to quantum behavior was understood to occur in 
the neighborjiood of a crossover temperature [3| . A few 
years later [4], Devoret et al. carried out another such 
experiment and again saw this evidence for MQT at the 
lowest temperatures attained. 

Beginning in 1985, Martinis, Devoret and Clarke [1, Q 
applied these ideas to the interpretation of experiments 
which added microwaves to the fixed current bias of an 
isolated Josephson junction. The thinking was that be- 
cause these experiments were carried out at low tempera- 
tures (18m_R' and 28mif ), only a quantum model need be 
considered. Biasing the junction very close to its critical 
current makes the well very shallow. Within the concep- 
tual framework of quantum phenomena, a well must con- 
tain discrete levels, and a shallow well will contain only a 
few of them. The action of properly tuned microwave ra- 
diation was then understood to populate the upper level 
from which escape would occur via macroscopic quantum 
tunneling of the virtual particle (i.e., the system would 
switch to the running state). The preferred frequencies 
for enhanced escape observed in the experiment were thus 



taken as signatures of the spacing of the discrete levels. 
The experimental data seemed to confirm this model Q 
based on the quantized one-dimesional Josephson wash- 
board potential. This quantum picture was embraced by 
the community of researchers interested in the possibil- 
ity of using Josephson based circuits as superconducting 
qubits. 

However, no classical considerations were made in [HI, 
so the argument as to an applicable model was in truth 
incomplete. In fact, as we show here, the classical RCSJ 
model, which is the foundation for the above-mentioned 
quantum phenomenology, with dc bias and applied mi- 
crowaves predicts the experimental observations to high 
accuracy. 

The phase dynamics of a Josephson junction which is 
biased with both dc and ac currents are governed by the 
simple classical expression [l| 



(/3 -I- q;(/3 -I- sin = 77 -I- ?7rf sin(il(iT) , 



(1) 



where the normalized dissipation is a = hu}o/2eRIc, 
with the Josephson zero-bias plasma frequency wq = 
y^2eIc/fiC. The junction RCSJ parameters are: critical 
current /c, capacitance C, resistance R. In this equation, 
the fixed and ac bias currents {ij and i]d) have been nor- 
malized to /c, the ac frequency has been normalized to 
Wq, and time is in units of uj^^ . Equation ([1]) can be 
viewed as describing the classical motion of a fictitious 
particle on a one dimensional surface 



U {(p) — —rjip — cos ip , 



(2) 



with the potential expressed in units of the Josephson 
coupling energy Ej = hlc/2e. This washboard potential 
has a repeated series of wells which become shallower as 
the bias 77 is increased. The natural frequency of small 
oscillations at the bottom of each well changes with bias 
according to Wp = a;o(l — 77^)4 . The position of the well 
minimum is tpmin = sin^^ 77 and the rim is located at 

In the absence of ac excitation (77^ = in Eq.((T])), 
the particle would sit at the bottom of the well. Suffi- 
cient additional thermal noise can activate escape over 



2 



1~\ ^ 0.9955 1 


\ Clp (a) 







11 = 0.9915 




(b) 



r]^0.9S95 




(c) 
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FIG. 1: Plots of maximum amplitude of the junction phase 
oscillation on the barrier side of the potential well as a func- 
tion of the normalized frequency of the applied microwave 
excitation for three different normalized bias currents i] (nu- 
merical simulation). In each case, the bias dependent nor- 
malized plasma frequency, determined from the expression 
flp = (1 ~ r]^)i , is shown by a dashed vertical line. The an- 
harmonic result fip, from Eq.Q is shown by a solid vertical 
line. The normalized microwave frequencies from the experi- 
ments described in are marked by inverted black triangles. 



the barrier AU = 2 
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tivation {rjd > 0) provides an alternate mechanism for 
escape, but the anharmonic shape of the well means that 
the ideal condition is not simply to set the normalized 
forcing frequency fid to {1 — rj^)* ; generally speaking a 
slightly lower value is optimum. It is also true that when 
the forcing frequency is properly tuned to a specific well, 
there will be some forcing amplitude above which escape 
will always occur. This means that just below such a 
drive threshold and when combined with even very small 
thermal noise, stochastic activation will occur, as is seen 
in all the experiments. 

A simple anharmonic analysis of the classical model 
can produce explicit results that can be compared di- 
rectly with experiments. This was first done in Refs. 
[^, where direct comparisons between classical theory 
and experiments showed that microwave-induced reso- 
nant switching can be understood classically. The key re- 
sults of the analysis come from introducing a monochro- 
matic ansatz if = ipo + asinQp^T (tpo is a constant and a 
is the oscillation amplitude) into Eq.(IT]) for a = ?7d = 0. 
This yields the effective relationships Jq (a) sin ipo — i] 
and 
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where J„ is the Bessel function of nth order and first kind. 
Since the experiments are conducted such that high es- 
cape probability exists near the sought-after resonance 



while vanishing escape probability should exist when the 
system is off-resonance, the strength of the microwave 
field should be such that the resulting oscillation ampli- 
tude reaches the inflection point of the effective potential; 
i.e., 
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Oscillating with this amplitude will result in a good sta- 
tistical probability that the junction switches to the non- 
zero voltage state within some allowed time. A useful 
explicit second order approximation in a for 77 — ?> 1 to 
the desired amplitude is given by — |(1 — 77). In- 
serting this result into Eq.Q provides an explicit rela- 
tionship between the anharmonic classical resonance and 
the bias current for direct comparisons with reported ex- 
periments in the limit of validity for Eq.Q. The pro- 
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TABLE I: Summary of experimental data for the peaks in the 
escape rates as reported in Martinis et al. [5]. 

tocol adopted in the experiments Q was to fix the 
microwave frequency and then observe the escape rate 
while scanning the bias current. Our attention here is on 
the results originally presented in most particularly 
the data in their Figs. 2 and 3. Two different physical 
samples were used - the parameters for Fig. 2 were given 
as /c = 30.572/iyl, C = 47.0pF, yielding a plasma fre- 
quency = 4.443 X 10^°s^^, whereas for Fig. 3 they 
were 1^ = 9.489^^ , C = 6.35pF, giving a plasma fre- 
quency ujQ — 6.734 X 10^°s~^. From blowups of these 
figures, it was possible to identify the peak positions, 
that is the bias current at maximum escape rate for each 
stated microwave frequency - these are given in the sec- 
ond and third columns of Table I. Normalized values fi^ 
and r\ are listed in the fourth and fifth columns. The 
most prominent peak in Fig. 2 and all four peaks in 
Fig. 3 are included. It should be noted that because 
the precise shapes of the potential wells are completely 
determined by the bias current 77, it is somewhat more 
physically meaningful, in contrast to the experimental 
protocol mentioned above, to fix the bias current and 
vary the excitation frequency. This approach reveals the 
optimum excitation frequency for any given well (bias). 
Such a strategy was followed in carrying out numerical 
solutions of Eq.([T]). In these simulations, the value of the 
dissipation constant a was chosen to be 0.00845, other- 
wise there were no adjustable parameters. For any se- 
lected bias 77 matching one of the experimental values, a 
series of runs was carried out for specific frequencies 17^. 
The excitation amplitude ry^ was ramped up from zero to 
its final value (typically in the range 0.0007 to 0.001) over 
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FIG. 2: A comparison of the classical and quantum models 
for the escape to a finite voltage state of a low temperature 
Josephson junction biased very close to its critical current and 
with applied microwaves. The 5 experimental data points are 
from Martinis et al. as is their predicted quantum result 
(dashed line) . The anharmonic analysis of the classical model 
yields the solid line. 



an interval of 600 plasma periods; this was done to avoid 
spurious start up effects. An additional 400 plasma peri- 
ods were then recorded. In each run the resulting data set 
from the numerical solution (^(r) was examined for the 
largest repetitive maximum angle reached. These maxi- 
mum angles are plotted in the three panels of Fig. 1. They 
correspond to the first, second and fifth rows of Table 
I. Computations corresponding to the other bias values 
yield the same general behavior. The obvious asymme- 
tries of the responses to ac excitation are a result of the 
anharmonicity of the wells. The inverted filled triangles 
mark the frequencies of the experimentally determined 
peaks in the escape rates, as reported in These fre- 
quencies coincide almost exactly with the location of the 
sudden onset of large induced phase oscillations in the 
numerical solutions of Eq.(IT]). These plots also reveal 
that good agreement should not be expected with the 
normalized linear plasma frequency Vlp because it is al- 
ways an over-estimate of the true optimum excitation for 
which a better expression is given by Eqs. ^ and 

Figure 2 presents a direct comparison of the classical 
and quantum models with respect to the experimental 
data reported in . The left-most four points come from 
the four peaks presented in Fig. 3(a) of Q; the fifth 
data point, on the far right, comes from the position of 
the right-most peak in Fig. 2(a) of Q. The dashed line 
labeled "Quantum" is obtained from the predicted energy 
level spacing Eq^i shown as a solid curve in Fig. 3(b) of 
01; the solid curve labeled "Classical" is generated from 
Eqs. dSl) and dH). 

As has been demonstrated, in the absence of noise, op- 
timum ac excitation will induce a maximum excursion 
in the phase of the Josephson junction. These largest 
excursions take the system closest to the critical value 
ificrit, at which escape would occur. Our basic hypoth- 
esis is simple: for a slightly subcritical excitation, ad- 
ditional noise will stochastically induce escapes, as was 
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FIG. 3: Switching distributions for swept bias i]. For these 
simulations, we chose ksT/Ej = 3.8418 x 10"^ and = 
0.282809 to match the 2 GHz frequency in Fig. 2 of [5|. Each 
distribution is generated from 100,000 switching events, and 
each simulation is initiated at 77 = 0. The applied sweep rate, 
microwave amplitude and damping were not reported in 
we have chosen a — 10~* and r]d = 0.0019. The applied sweep 
rate rj is reported on each panel. 



demonstrated in [9|. Therefore the experimental obser- 
vations of peaks in the escape rates at finite temperature 
are fully consistent with the classical model embodied in 
Eq. ^ . Agreement between predictions of the classical 
model with experimental values for the dc bias (77) and 
microwave frequency (fid) at these peaks is very accu- 
rate. Moreover, the classical model maintains its accu- 
racy over two distinct experimental samples having quite 
different critical currents and capacitances. One curious 
observation made in Q is the appearance of multiple res- 
onant escape peaks seen in their Figure 2. These were 
interpreted as evidence for the higher quantum transi- 
tions 1 — > 2, 2 ^ 3, and so on. However, no signatures 
of these were seen in the comparable escape probability 
plots shown in their Figure 3. We further notice that such 
a family of peaks was also not observed in [l^] . Our sim- 
ulations of the classical model show that closely spaced 
multiple resonant escape peaks can be observed under 
certain conditions due to the nonlinearity that provides 
the possibility for phase-locking of an oscillator at many 
harmonics (sub, super, and fractional). We mention in 
passing that this has also been shown ^ to explain the 
so-called multi-photon absorption peaks discussed in flo| . 
One could imagine that the multiple resonances seen in 
Figure 2 of Q could be due to some fractional resonance, 
but we see no evidence for this in the classical model for 
the given parameters. Instead, we have noticed very sim- 
ilar families of resonant peaks appearing when the bias 
current sweep rate is elevated. This important parameter 
is not given in Q , but typical values would correspond to 
a normalized sweep rate in the interval between ^ 10~^ 
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and ~ 10~^. It is in the upper end of this range (faster 
sweep rates) that we encounter the appearance of the 
multiple peaks similar to those observed in Q. Figure 3 
shows the appearance of such a switching peak family as 
a function of the bias sweep rate fj for a fixed microwave 
amplitude . All other parameters are chosen according 
to the experiments leading to figure 2 in Q. Thus we 
included in the right hand side of Eq.(IT]) a stochastic 
current noise term Q corresponding to the thermody- 
namic temperature of 2SmK\ each distribution in Fig. 3 
is the result of 100,000 switching events after initiating 
the system at = 0. 

We clearly observe in Fig. 3 the emergence of multiple 
peaks, and we further observe that the spacing between 
the peaks increases with the sweep rate. We have inves- 
tigated the details of the resonance near the bias values 
of the multiple peaks, and we have found only a single 
resonance, which is well described by Eqs. and 
and this is consistent with the fact that only a single 
resonant peak is found for small sweep rates. Thus, the 
multiplicity of apparent resonant peaks must be related 
to an interference between the microwave oscillation and 
the non-negligible variation in j] for high sweep rates. In 
this context it is worthwhile recalling that a fast bias 
current sweep rate has previously been shown to pro- 
duce multi-peak switching distributions in simulations of 
long Josephson junctions JTj, and it has been shown ex- 
perimentally that fast sweep rates can produce multiple 
peaks in the conditional switching distribution [T2| even 
without the application of microwaves. 

To conclude, we have established that the postulated 



quantum ingredients of discrete levels and macroscopic 
quantum tunneling are superfluous in explaining the 
overall results of experiments on resonant escape in cur- 
rent biased Josephson junctions at low temperature and 
with applied microwaves; the well-known classical model 
has been shown to yield excellent agreement with the re- 
ported experimental data. This agreement is consistent 
with the ones previously reported in Refs. [8, 9], and 
is quite similar to our recent examination of claims of 
quantum entanglement in_coupled Josephson qubits [31 , 
and to our earlier work pM^ showing that the classical 
model accounts for reports of, e_^., Rabi oscillations in 
large Josephson junction qubits jl5| . 

The often repeated claim Tb'] that the phase variable 
of Josephson junctions in these superconducting systems 
at these temperatures has been "unequivocally" demon- 
strated to be exclusively a quantum variable is not sup- 
ported by the present evidence. What is lacking are ex- 
perimental observations of phenomena that cannot be en- 
compassed within a classical picture and thus require an 
exclusively quantum model. The classical RCSJ model 
seems adequate in seeking an overall understanding of 
these types of experiment where microwaves are applied 
to Josephson systems. 
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